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A FLUID  MODEL  FOR  ESTIMATING  MINIMUM  SCALE  SIZES  IN 
IONOSPHERIC  PLASMA  CLOUD  STRIATIONS 


i 
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I.  INTRODUCTION 

Artificial  plasma  clouds  (i.e.,  nuclear  or  barium)  detonated  in 
the  ionosphere  typically  develop  visible  striations  over  periods  of  tens  of 
seconds  to  tens  of  minutes.  Propagation  studies  carried  out  simultaneously 
with  barium  cloud  experiments  during  project  STRESS  have  shown  that  line-of- 
sight  radio  communication  through  a striated  region  may  be  subject  to  10-30  db 
of  attenuation^".  Thus  local  enhancements  in  the  ionospheric  plasma  density 

can  have  a great  impact  on  propagation  when  small  scale  structuring  occurs. 

2 

The  resolution  of  the  optical  data  from  barium  cloud  experiments  and  nuclear 
detonations  is  insufficient  to  reveal  structures  (irregularities)  smaller 
than  100's  of  meters.  However,  recent  in  situ  rocket  measurements ^ suggest 
that  structure  down  to  meter  sizes  may  be  present.  We  wish  to  address 
theoretically  the  question  of  what  determines  the  smallest  structures  to  be 
found  in  a striated  plasma  cloud. 

The  mechanism  by  which  striations  are  produced  in  barium  clouds^-^ 
and  most  nuclear  clouds  has  been  established  to  be  the  gradient  drift 
instability.  This  is  an  electrostatic  fluid  process  analogous  to  the  Rayleigh- 
Taylor  instability.  The  results  of  the  model  are  in  need  of  correction  for 
kinetic  effects  below  scale  sizes  approximating  the  ion  gyro  radius.  For  a 
barium  ion  in  the  daytime  F region,  the  gyro  radius  is  approximately  7 meters. 
Note:  Manuscript  submitted  August  24,  1978. 
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Before  carrying  out  the  kinetic  corrections  (to  be  reported  elsewhere)  one 
should  ascertain  whether  or  not  the  fluid  model  tends  to  generate  structures 
on  such  a small  scale. 

II.  THE  ONE  LEVEL,  TWO  DIMENSIONAL  FLUID  MODEL 
For  barium  clouds  released  at  altitudes  of  approximately  200  km 
or  greater,  electron  and  ion  collision  frequencies  are  small  compared  to  their 

7 

gyrofrequencies.  As  a result  the  electrical  conductivity  of  the  plasma  is 


dominated  by  the  scalar  (Pedersen)  component, 

v. 

v _ ec  in 

z ~ n r 


where  n,  e,  c,  B,  v.  , and  Jb  are  respectively  the  ion  number  density, 
electronic  charge,  speed  of  light,  magnetic  field  strength,  ion-neutral 


(1) 


collision  frequency,  and  ion  gyrof requency . Assuming  all  parameters  in 
(1)  to  be  constant  except  n,  I obeys  a continuity  equation, 

|i  = _ y.(y  Z)  + K V2  I,  (2) 

where  V is  the  plasma  drift  velocity,  and  K is  a diffusion  constant  used 
to  model  the  effects  of  electronic  collisions.  The  value  of  K.  will  be 
estimated  later.  Let  us  adopt  Cartesian  coordinates  (x,y,z),  with  z in  the 
direction  of  the  (constant)  magnetic  field.  We  assume  that  all  variables 
depend  only  on  x and  y (striations  are  observed  to  be  mainly  magnetic  field 
aligned) . We  express  the  electric  field  as 

E (x,y)  = Eq  y - V«S  (x,y)  (3) 

where  Eq  is  the  constant  ambient  electric  field  and  <f>  is  the  potential  due 
to  the  presence  of  the  cloud.  It  is  convenient  to  use  a frame  moving  with 
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the  ambient  plasma  drift  velocity,  Eq  c/B  x.  In  this  frame,  we  have  a net 
plasma  drift  velocity 

V = - c/B  7/  x z (4) 

The  equation  can  be  viewed  as  a solution  of  the  electron  momentum  equation 
in  which  inertia  and  collision  terms  have  been  neglected.  The  collisional 
correction  has  been  included  in  an  approximate  way  in  the  diffusion  term  in 
(2).  The  set  of  fluid  equations  is  closed  in  this  simple  model  by  quasi- 
neutrality; i.e.,  the  constraint  that  the  electric  current  be  divergenceless: 

7*  (17/)  = E -7 l (5) 
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The  more  sophisticated  two  layer  model  allows  polarization  currents  generated 
by  the  cloud  to  flow  along  the  magnetic  field  direction  and  close  in  the  lower 
ionosphere,  where  cross-field  conductivities  are  high.  However,  the  lower 
layer  becomes  unimportant  when  the  cloud's  conductivity  is  high  enough  to 
let  current  loops  close  in  the  cloud  itself.  Some  of  the  barium  releases 
at  altitudes  of  approximately  200  km  are  known  to  have  integrated  conduc- 
tivities as  high  as  15  times  that  of  the  ionosphere.  For  such  large  clouds, 
the  one-level  model  should  be  adequate.  The  one-level  model  is  simple 
enough  to  allow  a convenient  scaling  law  to  emerge.  The  small  number  of 
variables  also  expedites  numerical  solution  of  the  equations  of  motion. 

III.  SCALING  THE  EQUATIONS  OF  MOTION 
Equations  (2),  (4),  and  (5)  can  be  put  into  dimensionless  form 
as  follows.  Let 

x = L x' 

- o - 

t = t t' 
o 

E*n' 

o 

V = V V', 

/ = LoEo/'  (6) 
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where  L is  an  arbitrary  scale,  length,  and 


V = cE  /B, 
o o 

t = L /V 

o o o 


(7) 


In  (6)  quantities  with  zero  subscript  are  dimensional  constants  and  all 
primed  variables  are  dimensionless.  Upon  expressing  (2),  (A),  and  (5)  in 
terms  of  the  dimensionless  primed  variables  and  dropping  primes  from  x',  t', 
I',  V’ , and  we  have 


~~  = -V*(VZ)  + k.'v2z 


-vrf  x 2 


where  K'  = 


V*(  EV<<)  = y • VZ 
K 


L V 
o o 


(8) 

(9) 

(10) 

(li) 


Note  that  K'  is  just  the  inverse  of  an  effective  Reynolds  number,  with  the 
diffusion  constant  K in  place  of  the  usual  kinematic  viscosity.  The  fact 
that  all  quantities  in  (8)-(10)  are  dimensionless  means  that  physical  systems 
of  different  sizes  will  evolve  in  the  same  way,  providing  initial  conditions 
and  boundary  conditions  are  analogous,  and  K'  has  the  same  value. 

Thus  the  model  (8)-(10)  can  be  used  to  answer  the  following  crucial 
question:  What  criterion  determines  whether  or  not  structures  in  a striated 
plasma  cloud  will  develop  into  structures  of  smaller  scale?  It  becomes  clear 
from  (8)  and  (11)  that  as  the  length  scale  becomes  smaller,  diffusion  becomes 
more  important.  For  sufficiently  small  scales,  we  expect  an  equilibrium  to 
be  reached  between  the  tendency  toward  finer  structuring  and  the  smoothing 
out  effect  of  diffusion.  Our  approach  to  finding  this  scale  will  be  to  carry 
out  a set  of  numerical  simulations  based  on  ( 2) — (5)  or  equivalently  (8)-(10) 
and  hopefully  to  identify  an  approximate  value  for  Kf  at  which  diffusion  is 


m . ■■■ 


just  able  to  halt  further  structuring  of  previously  formed  striations. 
Knowledge  of  a critical  value  of  K'  combined  with  estimates  for  the  drift 
speed  and  diffusion  coefficient  allow  calculation  of  a minimum  length  scale 
from  (11).  Drift  speeds  are  fairly  well  known  to  be  of  order  100  m/sec. 
However,  there  are  no  direct  measurements  of  the  diffusion  constant,  so  it 
must  be  calculated. 

IV.  CALCULATION  OF  THE  DIFFUSION  CONSTANT 
Neglecting  ion  and  electron  inertia,  the  momentum  equations  for 
electrons  and  ions  are 

kT. 

0 - n.V,  x Z - v,  V.  + v.  (V  - v, ) + - E - — - Vn/n  (12) 

i-i  in-i  le  -e  -l  m.  - m. 

i l 

kT 

0 = - 0 V xz-v  V +v,(V.  -V)-— E-  — ~ Vn/n  (13) 
e-e  en-e  ei  -l  -e  m - m 

e e 

For  the  purposes  of  this  derivation  we  have  adopted  a coordinate  system  at 

rest  with  the  neutral  atmosphere.  In  (12)  and  (13),  the  V's  are  fluid 
velocities,  with  subscripts  i and  e referring  to  ions  and  electrons. 

Similarly,  v stands  for  momentum  transfer  collision  frequency,  T for 
temperature,  and  m for  particle  mass.  The  gyrofrequencies  are 


“i  miC 

n = — (15) 

e m c 
e 

Note  that  we  have  retained  the  electron-ion  collision  frequency  in  (13) . 

2 -1 

For  conditions  typical  of  a 200  km  release,  ven  is  of  order  10  sec  , and 
3 -1 

is  of  order  10  sec  . This  is  because  the  ions,  although  less  numerous 
than  the  neutrals,  interact  with  electrons  via  a Coulomb  cross  section. 
However,  the  v term  in  (12)  is  typically  much  smaller  than  the  term 
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because  of  the  low  mass  of  the  electron.  Neglecting  v 


we  find  from  (12) 


and  (13) 


1 -1 

nV,  = M,  (r  nE 7n) 

-i  i Tn^  - dk 


1 ^ A 1 

nV  = M (-  — nE Vn)  + v .M  nV, 

-e  e m — m ei  e —I 

e e 


where 


V1  -f‘n  “l)/  < + 

V-Ql  Vin ' 

' P P ' 


V = V + V , 

e en  ei 


Let  us  invoke  quasineutrality  in  the  form 


7*(n7±)  = 7* (nVg) , 


and  take  the  divergence  of  (16)  and  (17).  We  assume  all  parameters  except 

n,  V, , V , and  E to  be  constant.  We  also  assume  all  vectors  have  only  x and 
-i  -e  - 

y components.  We  may  then  eliminate  7*nE  between  (16)  and  (17),  and  find, 


to  lowest  order  in  v /ft  , 

e e 


where 


-V* (n7  ) = - f-  7* (nE*z)  + KV2n, 

— X D “ 


v k(T  + T ) An 
K ’ “e  + »„  Uin/Cil 


1 
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We  use  the  convenient  approximations  , 

vei  . (34  + 4.18  log10(Te3/ne))  n^T,*372  (24) 

v * 1.8  x 105  p (25) 

en  rn  N 7 

vin  = 4 x 10  10  nn,  (26) 

where  pn  and  n^  are  the  neutral  pressure  and  number  density,  T£  is  in  ok  and 

all  other  quantities  are  in  cgs  units.  For  typical  daytime  conditions  at  200  km, 

we  have  T = 950°K,  T.  = 700°K,  p = 8.64  x 10  ^ erg  cm  3,  n = 7.0  x 10^  cm  3, 
e i n n 

= 9 x 10^  sec  3 , and  (singly  ionized  Ba)  = 36  s 3.  We  know  that  for 
typical  barium  releases^’ ^ 10^  5 n^  $ 10  2 cm  3.  Thus  K is  expected  to  lie 
between  0.6  and  6 m /sec.  When  the  calculation  is  repeated ifiing  data  appro- 
priate to  170  km  and  150  km,  K changes  by  less  than  a factor  of  2.  This  is 
because  is  dominated  by  ionic  collisions,  and  is  insensitive  to  moderate 
changes  in  the  neutral  background. 

V.  NUMERICAL  SIMULATION 

In  order  to  estimate  a critical  value  for  K',  we  have  solved 
numerically  equations  (2)-(5)  using  diffusion  constants  K determined  so  as 
to  give  one  of  four  desired  K'  values.  The  results  to  be  presented  here 
were  obtained  using  a mesh  of  162  x 42  grid  points  with  constant  grid  spacing 
in  both  directions  of  310  meters  and  doubly  periodic  boundary  conditions. 

Equation  (2)  was  advanced  in  time  using  a flux-corrected  leapfrog-trapezoid 

scheme,  which  is  basically  second  order  in  time,  fourth  order  in  space. 

Second  order  centered  differences  were  used  in  (4)  and  (5) . The  elliptic 

9,10 

potential  equation  (5)  was  solved  iteratively  using  the  Chebychev  method 

2 

The  initial  condition  (t=0)was  Z = 1 + 10  e -(x/8km)  (1  + e(x,y)),  where 

e was  generated  from  a k”^  power  spectrum  and  random  phases.  The  root  mean 
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J 


square  value  of  e was  0.03.  (The  dimensionality  of  I may  be  scaled  out  of 

(2)  and  (5)  with  no  consequence).  B is  taken  to  be  0.5  gauss  and  Vq  = 

cE  /B  = 100  m/sec. 
o 

We  are  not  concerned  here  with  details  of  relating  an  initial 
condition  to  a final  state  achieved  after  a long  period  of  development. 
Rather  we  wish  to  determine  whether  or  not  a given  configuration  tends 
toward  spontaneous  generation  of  structure  smaller  than  that  which  is 
present.  Thus  we  view  the  time  evolution  to  be  presented  here  as  a 
continuous  sequence  of  initial  conditions,  hopefully  at  some  point  descrip- 
tive of  structure  in  a striated  cloud.  We  seek  a value  of  the  Reynold's 
number  just  low  enough  that  a long  finger-like  striation  will  not  evolve 
into  new  and  smaller  scale  structures. 

Since  we  hope  to  identify  a critical  value  of  the  Reynold's 


number , 


V L 

R = -°-P 
Re  - K 


= 1/K' 


we  have  carried  out  numerical  simulations  so  as  to  keep  Rg  constant  through- 
out a given  computer  run.  Throughout  each  run  presented  here,  Vq  is  constant, 

and  L is  determined  self-consistently  from  the  conductivity  profile: 
o 

2 2 1/2 

Lo(t)  = {Z  (E(x,y,t)/E(VL(x,y,t))}i//  (28) 

where  l without  the  argument  (x,y,t)  refers  to  summation  over  all  grid  points. 
This  definition  of  Lq  is  somewhat  arbitrary,  but  it  does  have  the  desirable 
feature  of  being  sensitive  to  small  scale  structure.  Since  Lq  changes  in 
time,  K must  change  if  Rg  is  to  be  held  constant.  Thus  during  a computer 
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run,  we  evaluate  L at  each  timestep  and  determine  K from  (27).  It  may  be 
o 

more  conventional  to  keep  K fixed  and  let  Rg  vary  with  Lq,  but  out  approach 

allows  a more  direct  identification  of  a critical  R . 

e 

The  results  of  four  computer  runs  (Rg  = °°,  800,  400,  and  300) 

are  shown  at  five  times  (t  = 0,  240,  360,  480,  and  720  sec)  in  Figures  1-5. 

All  four  cases  show  the  breakup  of  the  initial  slab  geometry  into  structure 

elongated  in  the  x direction.  However,  only  the  top  two  cases  Rg  = <*>  and 

800  show  a clear  tendency  toward  further  structuring.  The  Rg  = 400  case 

appears  marginal,  and  the  Rg  = 300  case  is  clearly  held  back  by  diffusion. 

Therefore  we  tentatively  propose  that  the  critical  Reynold's  number  for  the 

initial  condition  presented  here  (E  /£  . = 11)  is  approximately  400.  We 

max  min 

intend  to  carry  out  several  more  simulations  using  finer  resolution  and 

different  initial  conditions  to  determine  if  R = 400  is  in  fact  universal 

e 

to  the  short  scale  limit  of  the  striation  process.  Using  Rg  = 400,  and 
Vq  = 100  m/sec.  in  (27),  we  find  the  following  range  for  Lq: 

For  n:  10^  - 10^  cm  ^ (29) 

2 

and  K:  0.6  - 6 m /sec. 

L : 2. 4 - 24  m 
o 

VI.  CONCLUDING  REMARKS 

The  one  level,  two  dimensional  striation  model  presented  here 
predicts,  for  the  particular  set  of  initial  conditions  considered,  minimum 
length  scales  of  2.4-24  meters  for  large  barium  clouds  near  200  km  altitude. 
These  scales  are  small  enough  that  kinetic  effects  may  be  important.  The 
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kinetic  corrections  have  not  yet  been  calculated,  and  will  be  reported 
elsewhere.  However,  the  tendency  of  the  fluid  mechanism  to  generate 
structures  this  small  is  in  apparent  agreement  with  recent  project  STRESS 
rocket  measurements^  and  with  the  attenuation  (probably  due  to  strong 
scattering  or  diffraction)  observed  in  propagation  experiments1.  Further 
work  needs  to  be  done  to  determine  the  sensitivity  of  the  results  to 
initial  conditions.  One  should  also  investigate  the  effects  of  ion  inertia 
(which  may  be  important  at  higher  altitudes) , and  electrostatic  coupling  to 
the  lower  ionosphere  in  the  case  of  small  clouds. 
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Figure  1 


Initial  conditions.  Contours  of  constant  I are  boundaries 
between  light  and  dark  areas.  A maximum  Z of  11.32  relative  to  the  back- 
ground occurs  in  the  light  vertical  band  containing  the  most  noticeable 
bump.  The  vertical  (y)  extent  of  each  rectangular  box  is  12.4  km;  the 
horizontal  (x)  extent  is  49.9  km.  The  finite  difference  grid  is  42  by  162 
points  in  the  y and  x directions  respectively.  Boundaries  are  doubly 
periodic.  The  Reynold's  numbers  are  (top  to  bottom)  R = °°,  800,  400,  300 


Figure  4 

£ contours  at  t = 480  seconds 
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6565  ARLINGTON  BLVD 
FALLS  CHURCH,  VA  22046 
OICY  ATTN  H.  BLANK 
OICY  ATTN  JOHN  SPOOR 
OICY  ATTN  C.  NAIL 
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COMSAT  LABORATORIES 
L I NTH  1 CUM  ROAD 
CLARKSBURG,  MD  20734 
01CY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

01CY  ATTN  D.  T.  FARLEY  JR 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 

01CY  ATTN  H.  LOGSTON 

01CY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 

01CY  ATTN  J.  ROBERTS 
01CY  ATTN  JAMES  MARSHALL 
01CY  ATTN  C.  W.  PRETTIE 

FORD  AEROSPACE  & COMMUNICATIONS  CORP 
3939  FABIAN  WAY 
PALO  ALTO,  CA  94303 

01CY  ATTN  J.  T.  MATTINGLEY 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.  O.  BOX  8555 
PHILADELPHIA,  PA  19101 

01CY  ATTN  M.  H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.  0.  BOX  1122 
SYRACUSE,  NY  13201 

01CY  ATTN  F.  RE I BERT 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
01CY  ATTN  DAS I AC 
01CY  ATTN  DON  CHANDLER 
01CY  ATTN  TOM  BARRETT 
01CY  ATTN  TIM  STEPHANS 
01CY  ATTN  WARREN  S.  KNAPP 
01CY  ATTN  WILLIAM  MCNAMARA 
01CY  ATTN  B.  GAMBILL 
01CY  ATTN  MACK  STANTON 

GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 

HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

01CY  ATTN  G.  MILLMAN 
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GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.  0.  BOX  6770 
SANTA  BARBARA,  CA  93111 

01CY  ATTN  JOHN  ISE  JR 
01CY  ATTN  JOEL  GARBARINO 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

CALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.  N.  DAVIS  (UNCL  ONLY) 

01CY  ATTN  NEAL  BROWN  (UNCL  ONLY) 

01CY  ATTN  TECHNICAL  LIBRARY 

GTE  SYLVAN I A,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  MARSHAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
URBANA,  IL  61803 

01CY  ATTN  K.  YEH 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
801  S.  WRIGHT  STREET 
URBANA,  IL  60680 

CALL  CORRES  ATTN  SECURITY  SUPERVISOR  FOR) 
01CY  ATTN  K.  YEH 


INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY -NAVY  DRIVE 
ARLINGTON,  VA  22202 

01CY  ATTN  J.  M.  AEIN 
01CY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFHARD 
01CY  ATTN  JOEL  BENGSTON 


HSS,  INC. 

2 ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01CY  ATTN  DONALD  HANSEN 

INTL  TEL  8 TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  HD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEMRA 
01CY  ATTN  JOHN  DASSOULAS 

LOCKHEED  MISSILES  6 SPACE  CO  INC 
P.  0.  BOX  504 
SUNNYVALE,  CA  94088 

01CY  ATTN  DEPT  60-12 
01CY  ATTN  D.  R.  CHURCHILL 

LOCKHEED  MISSILES  AND  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01CY  ATTN  MARTIN  WALT  DEPT  52-10 
01CY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
01CY  ATTN  W.  L.  IM-lOF  DEPT  52-12 

KAMAN  SCIENCES  CORP 
P.  0.  BOX  7463 
COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

LINKABIT  CORP 
10455  ROSELLE 
SAN  DIEGO,  CA  92121 

0 ICY  ATTN  IRWIN  JACOBS 

lOWELL  RSCH  FOUNDATION,  UNIVERSITY  OF 
450  AIKEN  STREET 
LOWELL,  MA  01854 

0 ICY  ATTN  K.  BIBL 


M.I.T.  LINCOLN  LABORATORY 


P.  0.  BOX 

73 

LEXINGTON, 

MA  02173 

01CY 

ATTN  DAVID  M.  TOWLE 

01CY 

ATTN  P.  WALDRON 

01CY 

ATTN  L.  LOUGHLIN 

01CY 

ATTN  D.  CLARK 

MARTIN  [MARIETTA  CORP 
ORLANDO  DIVISION 
P.  0.  BOX  5837 
ORLANOO,  FL  32805 

01CY  ATTN  R.  HEFFNER 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CY  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOULE 
01CY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
01CY  ATTN  R.  W.  HALPRIN 
01CY  ATTN  TECHNICAL  LIBRARY  SERVICES 


j 


MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 

01CY  ATTN  P.  FISCHER 
0 ICY  ATTN  W.  F.  CREVIER 
01CY  ATTN  STEVEN  L.  GUTSCHE 
01CY  ATTN  D.  SAPPENF IELD 
01CY  ATTN  R.  BOGUSCH 
01CY  ATTN  R.  HENDRICK 
01CY  ATTN  RALPH  KILB 
01CY  ATTN  DAVE  SOWLE 
01CY  ATTN  F.  FAJEN 
01CY  ATTN  M.  SCHEI8E 
01CY  ATTN  CONRAD  L.  LONGMIRE 
01CY  ATTN  WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 
P.  0.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
01CY  ATTN  C.  E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

01CY  ATTN  W.  HALL 
01CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD. 

SANTA  MONICA,  CA  90404 

01CY  ATTN  E.  C.  FIELD  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASSIFIED  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  3027 
BELLEVUE,  WA  98009 

01CY  ATTN  E.  J.  FREMOUW 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  1069 
BERKELEY,  CA  94701 

01CY  ATTN  A.  THOMPSON 
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R S D ASSOCIATES 

P.  0.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

01CY  ATTN  FORREST  GILMORE 
01CY  ATTN  BRYAN  GABBARD 
01CY  ATTN  WILLIAM  B.  WRIGHT  JR 
01CY  ATTN  ROBERT  F.  LELEVIER 
01CY  ATTN  WILLIAM  J.  KARZAS 
01CY  ATTN  H.  ORY 
01CY  ATTN  C.  MACDONALD 
01CY  ATTN  R.  TURCO 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

01CY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDROZIAN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

0 ICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.  0.  BOX  2351 
LA  JOLLA,  CA  92038 

01CY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  D.  SACHS 
01CY  ATTN  E.  A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 

SCIENCE  APPLICATIONS,  INC. 

HUNTSVILLE  DIVISION 
2109  W.  CLINTON  AVENUE 
SUITE  700 

HUNTSVILLE,  AL  55805 

01CY  ATTN  DALE  H.  DIVIS 

SCIENCE  APPLICATIONS,  INCORPORATED 
8400  WESTPARK  DRIVE 
MCLEAN,  VA  22101 

01CY  ATTN  J.  COCKAYNE 

SCIENCE  APPLICATIONS,  INC. 

80  MISSION  DRIVE 
PLEASANTON,  CA  94566 
01CY  ATTN  SZ 


SRI  INTERNATIONAL 
335  RAVENS WOOD  AVENUE 
MENLO  PARK,  CA  94025 


01CY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
01CY  ATTN  G.  SMITH 
01CY  ATTN  L.  L.  COBB 
01CY  ATTN  DAVID  A.  JOHNSON 
01CY  ATTN  WALTER  G.  CHESNUT 
01CY  ATTN  CHARLES  L.  RINO 
01CY  ATTN  WALTER  JAYE 
01CY  ATTN  M.  BARON 
01CY  ATTN  RAY  L.  LEADABRAND 
01CY  ATTN  G.  CARPENTER 
01CY  ATTN  G.  PRICE 
01CY  ATTN  J.  PETERSON 
01CY  ATTN  R.  HAKE,  JR. 

01CY  ATTN  V.  GONZALES 
01CY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.  P.  BOQUIST 

TRW  DEFENSE  & SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90278 

01CY  ATTN  R.  K.  PLEBUCH 
01CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 

VISIDYNE,  INC. 

19  THIRD  AVENUE 

NORTH  WEST  INDUSTRIAL  PARK 

BURLINGTON,  MA  01803 

01CY  ATTN  CHARLES  HUMPHREY 
01CY  ATTN  J.  W.  CARPENTER 


IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
UNCLASSIFIED  ONLY 


PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE 

ADVANCED  RESEARCH  PROJECTS  AGENCY  (ARPA) 

STRATEGIC  TECHNOLOGY  OFFICE 
ARLINGTON,  VIRGINIA 

CAPT.  DONALD  M.  LEVINE 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

DR.  P.  MANGE 
DR.  R.  MEIER 

DR.  E.  SZUSZCZEWICZ  - CODE  7127 
DR.  TIMOTHY  COFFEY  - CODE  6700 
DR.  S.  OSSAKOW  - CODE  6780 
DR.  J.  GOODMAN  - CODE  7560 


SCIENCE  APPLICATIONS,  INC. 
1250  PROSPECT  PLAZA 
LA  JOLLA,  CALIFORNIA  92037 

DR.  D.  A.  HAMLIN 
DR.  L.  LINSON 
DR.  D.  SACHS 


DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL  LABORATORY 
NOAA 

BOULDER,  COLORADO  80302 

DR.  A.  GLENN  JEAN 
DR.  G.  W.  ADAMS 
DR.  D.  N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 


A.  F.  GEOPHYSICS  LABORATORY 
L.  G.  HANSOM  FIELD 
BEDFORD,  MASS.  01730 

DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.  M.  FORBES 
DR.  T.  J.  KENESHEA 
DR.  J.  AARONS 
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OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VIRGINIA  22217 

OR.  H.  MULLANEY 


COMMANDER 

NAVAL  ELECTRONICS  LABORATORY  CENTER 
SAN  DIEGO,  CALIFORNIA  92152 

DR.  M.  BLEIWEI SS 
DR.  I.  ROTHMULLER 
DR.  V.  HILDEBRAND 
MR.  R.  ROSE 


U.  S.  ARMY  ABERDEEN  RESEARCH  AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MARYLAND 

DR.  J.  HEIMERL 


COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 

DR.  T.  CZUBA 


HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MASS.  02138 

DR.  M.  B.  MCELROY 
DR.  R.  LINDZEN 


PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PENNSYLVANIA  16802 


DR.  J. 

S.  NISBET 

DR.  P. 

R.  ROHRBAUGH 

DR.  D. 

E.  BARAN 

DR.  L. 

A.  CARPENTER 

DR.  M. 

LEE 

DR.  R. 

DIVANY 

DR.  P. 

BENNETT 

DR.  E. 

KLEVANS 

i 


UNIVERSITY  OF  CALIFORNIA,  LOS  ANGELES 

405  HILLGARD  AVENUE 

LOS  ANGELES,  CALIFORNIA  90024 


DR.  F.  V.  CORONITI 
DR.  C.  KENNEL 


UNIVERSITY  OF  CALIFORNIA,  BERKELEY 
BERKELEY,  CALIFORNIA  94720 

DR.  M.  HUDSON 


UTAH  STATE  UNIVERSITY 
4TH  N.  AND  8TH  STREETS 
LOGAN,  UTAH  84322 

DR.  P.  M.  BAN<S 
DR.  R.  HARRIS 
DR.  V.  PETERSON 
DR.  R.  MEGILL 
DR.  K.  BAKER 


CORNELL  UNIVERSITY 
ITHACA,  NEW  YORK  14850 

DR.  W.  E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 
DR.  E.  OTT 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MARYLAND  20771 

DR.  S.  CHANDRA 
DR.  K.  MAEDO 


PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NEW  JERSEY  08540 

DR.  F.  PERKINS 
DR.  E.  FRIEMAN 


INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VIRGINIA  22202 
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UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA.  15213 

OR.  N.  ZABUSKY 
OR.  M.  BIONDI 


